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The interplay between the finite volume and finite lattice spacing is investigated using lattice QCD simulations 
to compute the Landau gauge gluon propagator. Comparing several ensembles with different lattice spacings 
and physical volumes, we conclude that the dominant effects, in the infrared region, are associated with the use 
of a finite lattice spacing. The simulations show that decreasing the lattice spacing, while keeping the same 
physical volume, leads to an enhancement of the infrared gluon propagator. In this sense, the data from P = 5.7 
simulations, which uses an a ~ 0.18 fm, provides a lower bound for the infinite volume propagator. 



I. INTRODUCTION AND MOTIVATION 

In the past years the Landau gauge gluon and ghost propa- 
gators have been studied using lattice QCD methods in pure 
Yang-Mills theories - see, for example, Bll 4l5ll and references 
therein. The main goal is to be able to compute the infrared 
propagators having control on the finite size and finite volume 
effects. 

The lattice simulations were performed using large physi- 
cal volumes, let us say above (3 fm) 4 , both for the SU(2) and 
SU(3) gauge groups. The largest volumes simulated so far be- 
ing ~ (27 fm) 4 for the SU(2) gauge group [6] and ~ (17 fm) 4 
for SU(3) yj]. To achieve such large volumes, the simulations 
use a relatively large lattice spacing, ~ 0.2 fm to be compared 
with the typical non-perturbative scale in QCD of ~ 1 fm. It 
is well known that the infrared propagator decreases as the 
physical volume is increased and the main concern about the 
recent simulations has been the control of the finite volume 
effects. For such large volumes, which are much larger than 
typical nuclear sizes, finite volume effects should be negligi- 
ble or smaller than other possible systematics. Less studied is 
the interplay between finite lattice spacing and finite volume 
effects. Indeed, given the large lattice spacing used to sim- 
ulate such large physical volumes, the interplay between the 
two effects can introduce a bias on the propagators. 

The aim of the present work is to investigate the interplay 
between finite lattice spacing and finite volume effects for the 
pure SU(3) Yang-Mills theory and, in particular, for the gluon 
propagator. In principle our conclusions should also be valid 
for the SU(2) simulations. Having control of all possible ef- 
fects is crucial to compare lattice results with the outcome of 
other non-perturbative techniques. Furthermore, being able to 
remove the finite size and volume effects, the propagators can 
be used to model more reliably non-perturbative physics. 

In order to understand how the finite lattice spacing and 
the finite lattice volume change the propagator, we perform 
various simulations with several lattice spacings and physi- 
cal volumes and compare the propagators. The simulations 
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show, once more, that the infrared propagator decreases when 
the physical volume of the simulation is increased. However, 
they also show that the finite lattice spacing effects have a 
much larger impact on the infrared propagator. The propaga- 
tors are qualitatively unchanged but we observe that a large 
lattice spacing underestimates the gluon propagator in the in- 
frared region. 

The simulations also show that we have now good control 
on the Landau gauge gluon propagator for momenta above 
900 MeV. Furthermore, the extrapolations to the infinite vol- 
ume discussed here suggest that we are able to provide a gluon 
propagator free of finite size effects for momenta above 400 
MeV. 

The paper is organized as follows. In section HI1 we discuss 
the lattice setup and renormalization procedure. In sectionHHl 
profiting from having various ensembles with various lattice 
spacings and volumes, the extrapolation of the zero momen- 
tum gluon propagator to the infinite volume limit is investi- 
gated. In section [TV] the lattice spacing effects on the infrared 
gluon propagator are investigated. In section [V] we model the 
propagator in order to be able to provide an extrapolation to 
the infinite volume. Finally, in section[VT]we resume and con- 
clude. 



n. LATTICE SETUP AND RENORMALIZATION 
PROCEDURE 

In the current work we report on simulations for the pure 
gauge SU(3) Yang-Mills theory using the Wilson action at 
several p values and for different physical volumes. Table U 
lists the various ensembles used to compute the gluon propa- 
gator. 

Thegauge configurations were generated with the MILC 
code 111611 using a combined Monte Carlo sweep of seven over- 
relaxation updates with four heat bath updates. For each en- 
semble the autocorrelation time was measured for the plaque- 
tte. For all the ensembles, the autocorrelation time was below 
10 combined Monte Carlo sweep. The gluon propagator was 
measured by a separation of ten autocorrelation times, with a 
minimum separation of 100 combined sweeps. For thermal - 
ization, we disregarded the first 500 combined Monte Carlo 
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p 


a(fm) 


1 /a (GeV) 


L 


La (fm) 


#Conf 


5.7 


0.1838(11) 


1.0734(63) 


44 


8.087 


55 


5.7 


0.1838(11) 


1.0734(63) 


36 


6.617 


100 


5.7 


0.1838(11) 


1.0734(63) 


26 


4.780 


132 


5.7 


0.1838(11) 


1.0734(63) 


18 


3.308 


149 


6.0 


0.1016(25) 


1.943(47) 


80 


8.128 


55 


6.0 


0.1016(25) 


1.943(47) 


64 


6.502 


121 


6.0 


0.1016(25) 


1.943(47) 


48 


4.877 


104 


6.0 


0.1016(25) 


1.943(47) 


32 


3.251 


100 


6.2 


0.07261(85) 


2.718(32) 


80 


5.808 


70 


6.2 


0.07261(85) 


2.718(32) 


64 


4.646 


99 


6.2 


0.07261(85) 


2.718(32) 


48 


3.485 


87 


6.4 


0.05449(54) 


3.621(36) 


80 


4.359 


52 



TABLE I: Lattice setup. 



sweeps. For the conversion to physical units we use the lattice 
spacing as measured from the string tension lfl7ll . 

The gauge links were rotated to the minimal Landau gauge 
by minimizing the function 



F[g} = 



1 



YR&Tt\W(x)~\ 



(1) 



where V is the number of lattice points, N c i = 4 the number 
of space-time dimensions, N c = 3 the number of colors and 
Ufj(x) the gauge transformed link. For the minimization of 
F\g\ over the gauge orbits we used an overrelaxation algo- 
rithm H, except for the p = 5.7 and V = 18 4 , 26 4 , 36 4 and 
P = 6.0 for V = 32 4 ensembles, where the configurations were 
rotated to the Landau gauge using a Fourier accelerated steep- 
est descent method 01811 . The quality of the gauge fixing was 
monitored by looking at 0, the lattice version of d ■ A" aver- 
aged over all lattice points and all colors; see, for example, 
|[l|] for details and definitions. The minimization was stopped 
when 9 < 10" B . 

For the computation of the gluon propagator we follow the 
definitions given in fH. Further, in the following we use the 
tree level improved momentum definition 



1/j : 



2 fun 
— sin I — 



0,1, 



2 



(2) 



where a stands for the lattice spacing and the number of 
lattice points in direction /j. The statistical errors on the prop- 
agators were evaluated with the jackknife method. 

The propagators computed with the ensembles listed in Tab. 
U will be compared with the large volume simulations per- 
formed by the Berlin-Moscow- Adelaide group 0]. Table HT1 
summarizes the data generated by this collaboration. Note that 
both sets in Tab. |I]and|II]use the Wilson action. Note also, that 
the data in Tab. [TT] reaches a physical volume which is about 
twice the largest volume reported Tab. Q] Moreover, the data 
reported in [4] use the lattice spacing computed from ro, while 
we use the lattice spacing measured from the string tension. 



p 


a (fm) 


I /a (GeV) 


L 


La (fm) 


5.7 


0.1838(11) 


1.0734(63) 


64 


11.763 


5.7 


0.1838(11) 


1.0734(63) 


72 


13.234 


5.7 


0.1838(11) 


1.0734(63) 


80 


14.704 


5.7 


0.1838(11) 


1.0734(63) 


88 


16.174 


5.7 


0.1838(11) 


1.0734(63) 


96 


17.645 



TABLE II: Lattice setup for the data taken from reference 0]. Note 
that the Berlin-Moscow-Adelaide data was rescaled to be coherent 
with our lattice spacing definitions - see text for details. 



The values are not compatible within errors and, therefore, to 
compare the Berlin-Moscow-Adelaide results with ours, we 
have rescaled their data accordingly. 

The simulations using the ensembles listed in table HI and HI1 
are associated with different lattice spacings. In order to be 
able to compare the various sets, one needs to renormalize the 
propagator data. As described in the next section, the renor- 
malization is performed at the scale fj — 4 GeV and, therefore, 
any differences should show up in the infrared region. 

The gluon data shown below has momentum cuts to reduce 
the lattice artifacts. For momentum above 1 GeV only the 
momentum which verify the conical plus cylindrical cuts ^ 
are included in our analysis. For momentum below 1 GeV, we 
have considered all the lattice data. 



A. Renormalization 



In order to renormalize the lattice data the bare gluon prop- 
agator is fitted to the 1-loop inspired functional form 



D(q 2 ) 



In 



(*) 



(3) 



where y = 1 3 /22 is the anomalous gluon dimension for pure 
SU(3) Yang-Mills theory, in the momentum range [qminAmax\- 
The renormalized gluon propagator is given by 



D{q 2 )=Z R D Lat {q 2 )., 



(4) 



where Zr is the renormalization constant and Duit{q 2 ) the 
bare lattice gluon propagator. The renormalization constant 
Zg is chosen such that 



D(q 2 ) 



1 

7 



(5) 



and here we will use fj — 4 GeV. This choice for the renormal- 
ization scale will allow a direct comparison with the extrapo- 
lation to the infinite volume of D(0) performed in JH]. 

In what concerns the fitting range, for the highest momen- 
tum q max , we will use the highest momentum achieved in the 
simulation, namely q max — A3 GeV for the P = 5.7 ensem- 
bles, q max = 7.8 GeV for the p = 6.0, q max = 10.9 GeV for 
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(3 = 6.2 and q max = 14.5 GeV for p = 6.4 data. The smallest 
fitting momentum q m j„ is chosen as the smallest momentum 
where % 2 jd.o.f. « 1 or is closer to one. Given that the renor- 
malization scale used throughout this work is /j = 4 GeV, we 
also require q m i„ to be smaller than 4 GeV. The % 2 /d.o.f. as 



function of q m j n can be seen in Fig. Q] 

The fitting range for the various ensembles considered in 
the present work are 



L- 1 tfmin tfmax L tfmin tfmax L tfmin tfnuix L tfmin Qmax 

p = 5.7 44 1.43 4.29 p = 6.0 80 3.25 7.77 p = 6.2 80 2.95 10.87 p = 6.4 80 2.70 14.50 

36 2.23 4.29 64 3.24 7.77 64 2.52 10.87 

26 2.44 4.29 48 2.50 7.77 48 2.12 10.87 

18 1.47 4.29 32 2.81 7.77 



with all momenta given in GeV. Proceeding in the same way, 
the q m i n used to compute Zr for the Berlin-Moscow-Adelaide 
propagator data is 3.38 GeV for the 64 4 lattice data, 3.23 GeV 
for 72 4 , 2.76 GeV for 80 4 , 3.19 GeV for 88 4 and 1.08 GeV 
for 96 4 gluon propagator data. 

The values of Zr/o 2 are reported in Fig. [2] It follows that 
Zr/o 2 is volume independent, at least within the statistical 
precision of the simulations reported here. Fitting the plotted 
data to a constant gives 

% =0.6501 (65) GeV 2 (6) 
a 

with a % 2 /d.o.f. = 0.50. The errors were computed assuming 
Gaussian error propagation. 

Despite the volume independence of Zr/o 2 , in the follow- 
ing, to compute the renormalized gluon propagator, instead of 
taking the central value of ©, we will use the central value 
obtained from fitting directly the lattice data to ([3j. In this 
way, we hope to take into account possible systematics from 
using a limited statistics in each simulation. The renormalized 
gluon propagator for all ensembles can be seen in Fig. [3] In 
Fig. [4] we show the renormalized propagators, for a selected 
set of data, with a higher resolution in momentum scale. 

Figure [3] shows that for a given lattice spacing the infrared 
gluon propagator is reduced as one goes towards larger phys- 
ical volumes. This is well known and has been reported many 
times. From the data at p = 5.7, the value of D(Q) decreases 



by a factor of ~ 1.35 when going from the smallest to the 
largest volume. 

For larger momenta the lattice propagator is less dependent 
on the lattice volume and, as Fig. |4]shows, for momenta above 
~ 900 MeV the lattice gluon propagator seems to be indepen- 
dent of the lattice volume, in the sense that all data sets define 
a unique curve for D(q 2 ). In this sense, one can claim that 
above such momenta the lattice propagator is free of lattice 
artifacts. 

In what concerns the infrared region, no clear sign of 
a turnover on D{q 2 ) is observed when approaching q = 0. 
Therefore, one can claim that the lattice data points towards a 
finite and non-vanishing value for the zero momentum gluon 
propagator. 



III. ON THE INFINITE VOLUME LIMIT OF D(0) 

Let us now discuss the infinite volume limit for D(0). The 
various simulations performed using different P values allow 
for independent extrapolations. For the extrapolation to the 
infinite volume we consider either 

D(0) = £+A.(0) or D(0) = ^+D40). (7) 
The fits of the lattice data for the different sets give 



p 


Doo(0) GeV 2 


% 2 /d.o.f. 


Doe(O) GeV 2 


t/d.o.f. 


5.7 


9.26 ± 0.13 


0.8 


8.43 ± 0.61 


2.6 


5.7 [4] 


7.35 ± 0.30 


1.8 


6.1 ± 1.4 


1.3 


6.0 


10.157 ± 0.097 


0.3 


8.79 ± 0.64 


1.7 


6.2 


9.980 ± 0.055 


0.0 


9.72 ± 0.25 


0.1 




(a) P = 5.7 
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(c) p = 6.2 

FIG. 1: % 2 /d.o.f. for the fits 



where the first columns of values refers to the extrapolation 
assuming a 1 /V dependence and the second columns assumes 
a 1/L linear function. The fits using combined data sets, i.e. 
all data or all p = 5.7 data, give a too large % 2 /d.o.f., meaning 
that they are not described by any of the expressions in (0. 

The fits and the linear extrapolation assuming a 1 /L behav- 
ior can be seen in Fig. [5] Although the fits have acceptable 
% 2 /d.o.f, as seen in Fig. [5] the linear extrapolation in 1/L 
does not provide a coherent picture of all the data sets. The 
situation does not improve when we consider a linear depen- 
dence in 1 /V. Note that the linear extrapolation in 1 /V gives 
larger values for D(0), when compared to the 1 /L extrapola- 
tion. 

From the above results, one can claim a D(0) in the range 6 
-10 GeV 2 , in good agreement with the linear extrapolation 
performed by JH], which used improved actions and gauge fix- 
ing, and where it was estimated a L\o(0) = 7.95(13) GeV 2 . 
Recall that the extrapolation performed in [8] assumed a l/V 
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u 
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(d) p = 6.4 

Eq. <[3]l as a function of q m i„. 



linear dependence for the lattice data. 

The analysis discussed so far supports a finite and non- 
vanishing D(0). Indeed, no turnover of D(q 2 ) is observed 
when approaching the zero momentum limit. However, as 
described in ||9|] the scaling analysis of the Cucchieri-Mendes 
bounds [7] and the ratios defined in ifioll do not exclude com- 
pletely the possibility of having a vanishing zero momentum 
propagator. 



IV. FINITE LATTICE SPACING EFFECTS 

The renormalization of the gluon propagator is devised to 
remove lattice spacing effects. As can be observed in Fig. 
|U the various renormalized lattice gluon propagators agree, 
within errors, in the ultraviolet region but not necessarily at 
low momenta. 
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f, = 6.0 
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FIG. 2: Renormalization constant for the different simulations. The 
data is volume independent and a fit to a constant gives Zr/o 2 = 
0.6514(64) GeV 2 . The continuum line is the central value of the 
fitted data and the dashed lines are the one standard deviation values. 



The finite lattice spacing effects can be investigated com- 
paring the renormalized gluon propagator computed using the 
same physical volume but different p values. For the simu- 
lations reported in Tab. U the following four sets have close 
physical volumes: 



be larger than the corrections associated with the finite physi- 
cal volume. Indeed, if, within each set of volumes, D(0) seem 
to become compatible within two standard deviations, the dif- 
ference between propagators for the infrared region is clearly 
beyond two standard deviations. Furthermore, Fig [6] shows 
that the coarser lattice simulations, i.e. those with smaller (3 
value, underestimate D(q 2 ) in the infrared region. Note, how- 
ever, that qualitatively the propagator is unchanged. In this 
sense, the large volume simulations performed by the Berlin- 
Moscow-Adelaide group provide a lower bound for the con- 
tinuum infrared propagator. 

In what concerns the zero momentum gluon propagator, 
the previous analysis suggests that D(Q) computed at p = 5.7 
gives a lower bound to its continuum value. For complete- 
ness, below we list the zero momentum gluon propagator for 
the largest physical volume associated with each P value sim- 
ulation 



p 


La (fm) 


D(0) (GeV~ 2 ) 


5.7 


8.087 


9.42 ±0.28 


5.7 


16.178 


7.72 ±0.20 


6.0 


8.128 


10.50±0.50 


6.2 


5.808 


10.05 ±0.29 


6.4 


4.359 


10.06 ±0.36 



(i) (P 


= 5.1,L 


= 3.308 fm), 


(P 


= 6.0,L 


= 3.251 fm), 


(P 


= 6.2,L 


= 3.485 fm); 


(ii) (P 


= 5.7,L 


= 4.780 fm), 


(P 


= 6.0,L 


= 4.877 fm), 


(P 


= 6.2,L 


= 4.646 fm), 


(P 


= 6.4,L 


= 4.359 fm); 


(iii) (P 


= 5.7,L 


= 6.617 fm), 


(P 


= 6.0,L 


= 6.502 fm); 


(iiv) (P 


= 5.1, L 


= 8.089 fm), 


(P 


= 6.0,L 


= 8.128 fm). 



The physical volumes for the simulations considered in each 
data set do not match perfectly. However, the results summa- 
rized in Fig. |3] show a very smooth dependence of the gluon 
propagator with the lattice physical volume. Therefore, we ex- 
pect that the conclusions drawn for comparing the ensembles 
within each of the data sets illustrate how the lattice spacing 
changes D(q 2 ). 

In Fig. |6]the data for the various sets are reported and com- 
pared. From Fig. |6]one can see that for q above ~ 900 MeV 
the lattice data define a unique curve, i.e. the renormalization 
procedure removes all dependence on the ultraviolet cutoff a 
for the mid and high momentum regions. 

For the lower momenta Fig. |6]shows that D(q 2 ) is not in- 
dependent of a. Further, comparing Figs. [3] and [6] it follows 
that, in the infrared region, and for the lattice sizes L — 3 fm or 
above, the corrections due to the finite lattice spacing seem to 



V. MODELLING THE GLUON PROPAGATOR AND 
REMOVING THE LATTICE ARTIFACTS 

As described in the previous sections, the renormalized 
gluon propagator depends on the physical lattice volume, for 
the same lattice spacing, and on the lattice spacing, for the 
same physical volume. In what concerns the infrared region, 
it is observed that the propagator decreases slowly when the 
lattice physical volume is increased, while keeping the same 
lattice spacing. On the other hand, when the physical volume 
is keep constant, the infrared propagator is enhanced when the 
lattice spacing decreases. According to the results described 
above, the dominant effect is associated with the lattice spac- 
ing rather than with the volume dependence. 

In the present work we aim to provide a gluon propaga- 
tor which is free of finite size effects. Such a goal can be 
achieved, at least partially, through modeling the propagator. 
Of the possible functional forms, we chose to fit the lattice 
data to 

. q 2 +M 2 

D(q 2 )=Z , , ' t, (8) 

KH ' q 4 +M 2 q 2 +M* ' 

where Z is dimensionless and Mi , Mi and M3 have dimen- 
sions of mass. The above expression is the tree level predic- 
tion of the so-called refined Gribov-Zwanziger action supplied 
by an extra fitting normalization parameter Z. For Z = 1 , as 
shown in lfl9ll . Eq. ([8]l describes the lattice gluon propagator 
in the deep infrared and up to momenta ~ 1 .5 GeV. Although 
in the framework of the refined Gribov-Zwanziger action the 
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Renormalized Gluon Propagator - (i = 4 GeV 

P = 5.7 data 



Renormalized Gluon Propagator - (i = 4 GeV 

P = 6.0 data 
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La = 
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La = 


4.780 fm 






La = 


6.617 fm 




A 


La = 


8.087 fm 




o 


La = 


1 1.763 fm 






La = 


13.234 fm 






La = 


14.704 fm 






La = 


16.174 fm 






La = 


17.645 fm 






Renormalized Gluon Propagator - (i = 4 GeV 

p = 6.2 data 




• 


La = 


3.251 fm 






La = 


4.877 fm 






La = 


6.502 fm 




A 


La = 


8.128 fm 





1 1.5 
q [GeV] 



Renormalized Gluon Propagator - (i = 4 GeV 

P = 6.4 data 




FIG. 3: Renormalized gluon propagator for /j = 4 GeV for all lattice simulations. 



different mass parameters have a precise meaning and are as- 
sociated with various condensates and the Gribov mass pa- 
rameter, here we do not explore such a connection. Besides 
the theoretical motivation coming from the refined Gribov- 
Zwanziger action, © can be viewed as a Pade approximation 
to the lattice data. 

The functional form ([8]) does not include any log terms gen- 
erated by the loop corrections. Therefore, we do not expect 
that it describes the full set of momenta. Being so, we have to 



define a fitting range, which includes the infrared region and 
goes beyond the 1 GeV, where all data sets define a unique 
curve, and where ([8j is able to reproduce properly the lattice 
propagator. Although in our analysis we have considered var- 
ious fittings ranges, in the following we will report only the 
results for a fitting range starting at q = and going up to 
<lmax = 4 GeV. The results for the other fitting ranges being 
similar. The outcome of the fits being 
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p 


L 


Z 


M\ 


M\ 


M\ 


X 2 /d.o.f 


5.7 


18 


0.8290(91) 


3.94(18) 


0.583(41) 


0.3224(89) 


1.86 




26 


0.8156(74) 


4.16(15) 


0.603(32) 


0.3533(78) 


2.32 




36 


0.8254(62) 


3.98(11) 


0.557(23) 


0.3523(58) 


1.96 




44 


0.8231(71) 


4.08(12) 


0.583(25) 


0.3569(64) 


1.65 


6.0 


32 


0.811(11) 


4.45(22) 


0.694(48) 


0.3243(94) 


1.15 




48 


0.7996(99) 


4.51(19) 


0.669(37) 


0.3384(83) 


1.81 




64 


0.8209(64) 


4.21(11) 


0.604(22) 


0.3362(50) 


1.42 




80 


0.8196(59) 


4.22(10) 


0.631(22) 


0.3177(41) 


0.81 


6.2 


48 


0.806(14) 


4.32(25) 


0.623(49) 


0.342(11) 


1.25 




64 


0.8135(75) 


4.36(14) 


0.653(28) 


0.3575(68) 


0.79 




80 


0.8169(96) 


4.38(17) 


0.654(34) 


0.3635(82) 


1.15 



where all the mass parameters are given in powers of GeV. 
For each P set we combine all volumes and perform a linear 
extrapolation in l/L to the infinite volume of each parameter 
independently. The % 2 /d.o.f for the Z, M 2 , M 2 , M\ asso- 
ciated with the extrapolation being, respectively, 0.78, 0.65, 
0.64, 1.01 for the p = 5.7 ensembles, 1.23, 0.56, 1.05, 4.61 
for the P = 6.0 ensembles and 0.00, 0.00, 0.07, 0.04 for the 
P = 6.2 ensembles. The extrapolated parameters being 



P = 

Z 

M 2 

M\ 

M\ 

D(0) 



5.7 

0.821(10) 

4.09(17) 

0.558(36) 

0.380(11) 

8.84(45) 



6.0 

0.830(13) 

4.01(16) 

0.565(46) 



6.2 

0.83333(17) 

4.473(21) 

0.704(29) 

0.3959(54) 

9.42(14) 



where D(0) is given in GeV -2 . The extrapolated propagators 
together with the largest volume simulations results for each 
P are shown in Fig. Q Note that due to the bad extrapolation 
for M 4 when using the p = 6.0 ensembles, Fig. [7] does not 
include any information coming from these data. 

From Fig. [7] comparing the extrapolated p = 5.7 propaga- 
tor using volumes up to (8. 1 fm) 4 with the largest volume sim- 
ulated by the Berlin-Moscow-Adelaide group, it follows that 
the two propagators are compatible within two standard devia- 
tions. The extrapolated propagator being above the (17.1 fm) 4 
lattice simulation results. This supports, once more, the pre- 
vious observation that p = 5.7 results provides a lower bound 
for the infinite volume propagator. Furthermore, the extrapo- 
lated propagator from the P = 6.2 data is above the extrapo- 
lated p — 5.7 propagator. Note, however, that the extrapolated 
D(q 2 ) are compatible within two standard deviations. 

In what concerns the infinite volume limit for D(0), the ex- 
trapolations give a value of 8.84(45) GeV 2 from the p = 5.7 
data and 9.42(14) GeV 2 from the p = 6.2 data. These values 
are slightly above the infinite volume limit Doo(0) = 7.95(13) 
GeV 2 estimated in OJ. 



VI. CONCLUSIONS 

In this work, the effects of using finite volume and finite 
lattice spacing to compute the Landau gauge gluon propaga- 
tor in lattice QCD simulations are investigated. The propaga- 
tor was computed for various physical volumes, all above (3 
fm) 4 , and various lattice spacings ranging from 0.18 fm down 
to 0.054 fm. Our simulations confirm, once more, that the 
infrared propagator decreases as the lattice volume increases. 
Furthermore, the comparison between the data generated at 
similar physical volumes but different lattice spacings, show 
that, in what concerns the infrared momenta, the finite lattice 
spacing effects are larger than the finite volume effects. Our 
analysis shows that decreasing the lattice spacing leads to an 
increase in D(q 2 ) in the low momenta region. In this sense, the 
data coming from the large volume simulations by the Berlin- 
Moscow-Adelaide group, which have an a ps 0.18 fm, pro- 
vides a lower bound for the continuum propagator. In princi- 
ple, the same behavior is expected for the large volume simu- 
lations performed for the SU(2) gauge group where a « 0.22 
fm. This observation is also supported by extrapolations of 
the lattice data to the infinite volume limit. 

Besides the finite lattice and finite volume effects, we also 
investigate the zero momentum gluon propagator. The various 
extrapolations performed here point towards a finite value for 

Z) (°)L=4GeV around9GeV_2 - 
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FIG. 4: Zoom of the renormalized gluon propagator for /j = 4 GeV. 
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FIG. 6: Comparing the renormalized gluon propagator at fj = 4 GeV for various lattice spacings and similar physical volumes. 
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FIG. 7: Extrapolation of the gluon propagator to the infinite volume 
limit for the (3 = 5.7 and (3 = 6.2 ensembles. The full line repre- 
sents the D(q 2 ) computed with {8]l and using the central values of 
the extrapolated parameters. The dashed lines are the one standard 
deviations for D(q 2 ) computed assuming gaussian error propagation. 



